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Abstract: The sensitivity of the atmospheric neutrino experiments to the non-standard 
flavor-dependent interaction in neutrino propagation is studied under the assumption that 
the only nonvanishing components of the non-standard matter effect are the electron and 
tau neutrino components e ee , e eT , e TT and that the tau-tau component satisfies the con¬ 
straint e TT = |e eT | 2 /(l + e ee ) which is suggested from the high energy behavior for atmo¬ 
spheric neutrino data. It is shown that the Superkamiokande (SK) data for 4438 days 
constrains |tan/3| = \e eT /(l + e ee )| < 0.8 at 2.5 a (98.8%) CL whereas the future Hyper- 
kamiokande experiment for the same period of time as SK will constrain as | tan /3| < 0.3 
at 2.5frCL from the energy rate analysis and the energy spectrum analysis will give even 
tighter bounds on e ee and \e eT \. 
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1 introduction 

From the experiments with solar, atmospheric, reactor and accelerator neutrinos it is now 
established that neutrinos have masses and mixings [1]. Neutrino oscillations in the stan¬ 
dard three-flavor scheme are described by three mixing angles, 6 12 , 0 13 , 6 23 , one CP phase 
<5, and two independent mass-squared differences, and Am^. The sets of the param¬ 

eters (A m, 2 i, $ 12 ) and (lAm^l, $ 23 ) were determined by the solar neutrino experiments 
and the KamLAND experiment, and by atmospheric and long baseline neutrino experi¬ 
ments, respectively. $13 was determined by the reactor experiments and the long baseline 
experiments [1]. The only oscillation parameters which are still undetermined are the value 
of the CP phase 6 and the sign of Am^ (the mass hierarchy). In the future neutrino long- 
baseline experiments with intense neutrino beams the sign of Am^ and 5 are expected to 
be determined [2, 3]. As in the case of B factories [4, 5], such high precision measurements 
will enable us to search for deviation from the standard three-flavor oscillations (see, e.g., 
Ref. [ 6 ]). Among such possibilities, in this paper, we will discuss the effective non-standard 
neutral current flavor-dependent neutrino interaction with matter [7-9] given by 

^eJ 1 = -2V2e f Q P p G F {V altl P L vp) (f^Pf), (1.1) 

where / and f stand for fermions (the only relevant ones are electrons, u and d quarks), 
Gp is the Fermi coupling constant, and P stands for a projection operator that is either 
Pl = (1 — 7s)/2 or Pn = (1 + 75)/2. If the interaction (1.1) exists, then the standard 
matter effect [7, 10] is modified. We will discuss atmospheric neutrinos which go through 
the Earth, so we make an approximation that the number densities of electrons ( N e ), 
protons, and neutrons are equal . 1 Defining e a p = + 3e“^ + 3e(^, the hermitian 

1 This assumption is not valid in other environments, e.g., in the Sun. 
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3x3 matrix of the matter potential becomes 


A = A 


1 "E ^ee 

tfie 

^re ^r/x 



( 1 . 2 ) 


where A = \f2GpN e stands for the matter effect due to the charged current interaction in 
the standard case. With this matter potential, the Dirac equation for neutrinos in matter 
becomes 


. d 


( v e (x) ^ 

= [U diag (0, A£ 2 1 ,A£ 3I ) U~ l + A] 

( "e(x) \ 

\ v t{x) j 


\ M x ) j 


where U is the leptonic mixing matrix defined by 


C 12 C 13 

U = I -S 12 C 23 - ci 2 S2?,size l& 


S 12 C 13 
C 12 C 23 — Sl2S23Sl3e 


Si3e 


-iS 


,i5 


S23C13 


•S12S23 — Ci2C23Sl3e 


i5 


iS 


-C12S23 — ■Sl2C23'Sl3e C 23 C 13 


(1.3) 


(1.4) 


and AEjk = Am 2 k /2E = (: m 2 — m 2 k )/2E, Cjk = cos 9jk, Sjk = sin 0j k . 

Constraints on e a p from various neutrino experiments have been discussed in Refs. [11- 
19]. Since the coefficients e a p in Eq. (1.2) are given by e Q( g ~ e o /3 + 36“^ + 3e^, considering 
the constraints by Refs. [11-19], we have the following limits [20] at 90%CL: 


^ |e ee | < 4 x 10° |e eM | < 3 x 10 1 |e eT | < 3 x 10° \ 

|e MM | < 7 x 10- 2 |e MT | < 3 x KT 1 . (1.5) 

y |e Tr | < 2 x 10 1 J 


From Eq. (1.5) we observe that the bounds on e ee , e er and e TT are much weaker than those 
on e afl (a = e,/x,r). 

On the other hand, the non-standard interaction (NSI) with components e Q( g (a, (3 = 
e, r) must be consistent with the high-energy atmospheric neutrino data. It was pointed 
out in Refs. [21, 22] that the relation 


Cer| — 6 tt (1 ~b £ee) i 


( 1 . 6 ) 


should hold for the matter potential (1.2) to be consistent with the high-energy atmospheric 
neutrino data, which suggest the behavior of the disappearance oscillation probability 

/Am 2 r\ 1 

1 - P(^ u, t ) ~ sin 2 2 d atm sin 2 ( ^ j oc — 2 , (1.7) 

where sin 2 20 atm and Am 2 tm are the oscillation parameters in the two-flavor formalism. In 
Ref. [23] it was shown that, in the high-energy behavior of the disappearance oscillation 
probability 

1 - P(. u n -A v^) c 0 + + O ^^ 2 ^ i (1-8) 
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in the presence of the matter potential ( 1 . 2 ), |co| <C 1 and |ci| <C 1 imply e e/Jj — ~ 

C^t — 0 and e rr — |cer| / (1 T Cee)* 

Taking into account the various constraints described above, in the present paper we 
take the ansatz 


A = A 


(l 


\ 


+ 


c er 

0 


(1.9) 


V(1 + €ee) 


and analyze the sensitivity to the parameters e a p (a, (3 = e, r) of the atmospheric neutrino 
experiment at Superkamiokande and the future Hyperkamiokande (HK) facility [24]. 

The constraints on e ee and e eT from the atmospheric neutrino has been discussed in 
Refs [25-27] with the ansatz different from ours. 

The effect of the non-standard interaction in propagation for solar neutrinos has also 
been discussed in Refs. [15, 18, 19, 28, 29], and Refs. [19] and [29] give a constraint —0.06 < 
sin 6*23 < 0.41 (at 90%CL) and |ef^| < 0.4 (at Ay 2 = 4 for 2 d.o.f.), respectively. 

The sensitivity of the ongoing long-baseline experiments to the non-standard interac¬ 
tion in propagation was studied for MINOS in Refs. [30-34], and for OPERA in Refs. [35, 
36]. As for the future long-baseline experiments, the sensitivity of the reactor and super¬ 
beam experiments was discussed in Ref. [37], that of the T2KK experiment was studied in 
Refs. [23, 38], and that of neutrino factories [ 6 ] was discussed by many authors [39-46]. 

The paper is organized as follows. In sect. 2, we analyze the SK atmospheric neutrino 
data and give the constraints on the parameters e a p (a, (3 = e, r) from the SK atmospheric 
neutrino data. In sect. 3, we discuss the sensitivity to e Q( g (a, (3 = e,r). of the future 
Hyperkamiokande atmospheric neutrino experiment In sect. 4, we draw our conclusions. 


2 The constraint of the Superkamiokande atmospheric neutrino experi¬ 
ment on e ee and |e er | 

In this section we discuss the constraint of the SK atmospheric neutrino experiment on the 
non-standard interaction in propagation with the ansatz (1.9). The independent degrees of 
freedom in addition to those in the standard oscillation scenario are e ee , |e er | and arg(e er ). 

The SK atmospheric neutrino data we analyze here is those in Ref. [47] for 4438 days. 
In Ref. [47], the contained events, the partially contained events and the upward going /j, 
events are divided into a few categories. Since we have been unable to reproduce all their 
results of the Monte Carlo simulation, we have combined the two sub-GeV ^-like data set 
in one, the two multi-GeV e-like in one, the two partially contained event data set and the 
multi-GeV /i-like in one, and the three upward going // in one. Ref. [47] gives information 
on the ten zenith angle bins, while that on the energy bins is not given, so we perform 
analysis with the ten zenith angle bins and one energy bin, i.e., we perform the rate analysis 
as far as the energy is concerned. 

The analysis was performed with the codes which were used in Refs. [48-50]. y 2 is 
defined as 

X (Xsub—GeV Xmulti—GeV T ^upward) ■ (^-1) 

023 ,|Am§ 2 |,< 5 ,arg(e er ) 
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In eq. (2.1) y 2 for the sub-GeV, multi-GeV, and upward going /i events are defined by 
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( 2 . 2 ) 


I'e) 


(2.3) 


Xupward = min AAA[p) "#*) + a u N ji u P ~+ v p) 

^ j = 1 ^ 


The summation on j runs over the ten zenith angle bins for each x 2 > n j{ a ) (a=s, m, 
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u; a = e,fi) stands for the neutrino and antineutrino data of the numbers of the sub- 
GeV, multi-GeV, and upward going /i events, Nj(u a —» up) (N^(u a — >• up)) stands for the 
theoretical prediction for the number of i'g-like events {Ip = e, p) which is produced from 
up (up) that originates from u a (u a ) through the oscillation process u a — >• up (u a — >• up), 
and it is expressed as the product of the oscillation probability P(u a —>• up) ( P(u a —>• up)), 
the flux F(u a ) ( F(u a )), the cross section, the number of the targets and the detection 
efficiency. a a (a = s, m, u) stands for the uncertainty in the overall flux normalization for 
the sub-GeV, multi-GeV, and upward going p, events, (3 a i (f5 a 2 ) stands for the uncertainty 
in the relative normalization between u e - u /t flux (u - u flux) for the sub-GeV (a = s) 
and multi-GeV (a = m) events, respectively. It is understood that x 2 is minimized with 
respect to a s , (3 s k (k = 1,2), a m , [3 m k (k = 1,2), a u . We have put the systematic errors 

CTp s \ — & pm 1 — 0.03, £7^2 — &Pm2 — 0.05, (Tq, — 0.2 (2.4) 

and we have assumed that a s and a m for the contained events are free parameters as in 
Ref. [51]. We have omitted the other uncertainties, such as the E u spectral index, the 
relative normalization between PC and FC and up-down correlation, etc., for simplicity. 
In Eq. (2.1) the sum of each x 2 is optimized with respect the mixing angle 623, the mass 
squared difference |AmJ 2 |, the Dirac CP phase 6 and the phase arg(e er ) of the parameter 
e er . The other oscillation parameters give little effect on y 2 , so we have fixed them as 
sin 2 26*12 = 0.86, sin 2 2 #i 3 = 0.1 and Am 21 = 7.6 x 10“ 5 eV 2 . 


<U 


OJ 



Figure 1 . The allowed region in the (e ee , |e eT |) plane from the SK atmospheric neutrino data for 
a normal mass hierarchy (left panel) and for an inverted mass hierarchy (right panel). 

The result for the Superkamiokande data for 4438 days is given in Fig. 1. The best-fit 
point for the normal (inverted) hierarchy is (e ee , |e er |) = (-1.0, 0.0) ((3.0, 1.7)) and the 
value of x 2 a t this point is 79.0 (78.6) for 50 degrees of freedom, and goodness of fit is 2.8 
(2.7) crCL, respectively. The best-fit point is different from the standard case (e ee , |e er |) = 
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(0, 0), and this may be not only because we have been unable to reproduce the Monte Carlo 
simulation by the Superkamiokande group, but also because we use only the information on 
the energy rate and the sensitivity to NSI is lost due to the destructive phenomena between 
the lower and higher energy bins (See the discussions in subsect. 3.1). The difference of the 
value of X' 2 for the standard case and that for the best-fit point for the normal (inverted) 
hierarchy is Ay 2 = 2.7 (3.4) for 2 degrees of freedom (1.1 erCL (1.3 <tCL)), respectively, 
and the standard case is certainly acceptable for the both mass hierarchies in our analysis. 
From the Fig. 1 we can read off the allowed region for | tan/3| = |e eT |/|l + e ee \, and we 
conclude that the allowed region for | tan (3 | is approximately 

| tan J eT I —t < 0.8 at 2.5crCL. 

11 T Gel 

3 Sensitivity of the Hyperkamiokande atmospheric neutrino experiment 

to e ee and \e eT \ 

In this section we discuss the potential sensitivity of HK to e ee and |e er |. Here we assume 
for simplicity that the the Hyperkamiokande detector has the same detection efficiencies 
as those of SK, and that the fiducial volume of HK is twenty times as large as that of 
SK. Since HK is a future experiment, the simulated numbers of events are used as “the 
experimental data”, and we vary e ee and e er as well as the standard oscillation parameters 
trying to fit to “the experimental data”. Here we perform an analysis on the assumption 
that we know the mass hierarchy, because some hint on the mass hierarchy is expected 
to be available at some confidence level by the time HK will accumulate the atmospheric 
neutrino data for twenty years. 

Since “the experimental data” are the simulated numbers of events, we can perform a 
energy spectrum analysis, assuming that the detection efficiency etc. are all equal among 
neutrinos and antineutrinos. Before we study the sensitivity to NSI, as a benchmark of our 
analysis, we have investigated the significance of the wrong mass hierarchy with our code, 
assuming the standard oscillation scenario and using different numbers of the energy bins. 
By comparing our result with the one in Ref. [24], we have found that our analysis on the 
mass hierarchy gives a result similar to that in Ref. [24], when we work with two energy 
bins in the contained events (the sub-GeV and multi-GeV events) and the systematic errors 
which are slightly different from those in Ref. [24]. We have checked that the sensitivity to 
NSI is not affected significantly by changing the systematic errors. As for the upward going 
// events, since our ansatz (1.9) is taken in such a way that the oscillation probability with 
e a/3 (a, /? = e, r) approaches to the one with the standard scenario in the high energy limit, 
the upward going // events are expected to give a small contribution to the significance of 
NSI. So in the case of the energy spectrum analysis we will work with two energy bins in 
the contained events and a single energy bin in the upward going [i events. 

3.1 The case with the standard oscillation scenario 

First of all, let us discuss the case where “the experimental data” is the one obtained with 
the standard oscillation scenario. The values of the oscillation parameters which are used 
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to obtain “the experimental data” are the following 2 : 

A 777,31 = 2.5 x 10~ 3 * eV 2 , sin 2 023 = 0.5 ,5 = 0 , 
sin 2 2012 = 0.86, sin 2 20 i3 = 0.1, Amiji = 7.6 x 10 _ 5 eV 2 . (3.1) 

As in the case of the analysis of the SK data, we vary the oscillation parameters 023, | Am 32 |, 
S and arg(e er ) while fixing the other oscillation parameters sin 2 2012 = 0.86, sin 2 20 i 3 = 0.1 
and Am 2 ! = 7.6 x 10~ 5 eV 2 . 

In the energy rate analysis, y 2 is the same as (2.1) where the numbers of events are 
calculated with the standard oscillation scenario with djk, A?ri 2 fc and 5 given by Eq. (3.1), 
and we have assumed that all the systematic errors except crp m 2 are the same as those in 
Eq. (2.4) in the analysis of SK data. crg m 2 = 0.16, which is the uncertainty in the relative 
normalization between the v - v flux, was chosen because this value was used in the energy 
spectrum analysis on the significance of the wrong mass hierarchy to give the result close 
to that in Ref. [24]. (See the discussions below.) 

In the spectrum analysis, on the other hand, y 2 u b_GeV an d Xmuiti-GeV are re pl ace d by 

2 To distinguish the oscillation parameters for the “the experimental data” (J — 1, ■ • • , 10; A = 

L,H;a = s,m;l = e,/r), etc.) and those for the numbers of events (N%j(v a —¥ vp) (j = 1, • • • , 10; A = 

L, H;a = s,m\a, (5 = e, fi), etc.) for fitting, the parameters with a bar denote those for “the experimental 

data”, whereas those without a bar denote the parameters for the numbers of events for fitting. 
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(3.3) 


In Eq. (3.3) we have introduced the relative normalization, which in general depends on 
the flavor and the energy of the events, between the upward and downward going bins: 


7ai,2 



7ai, 2 (j < job A = L.H) 
-7A1,2 (j > ith ; a4 = L,H) 

71,2 (j < jth) 

-71,2 (j > jth), 


and jth = 3 is the index which separates the upward and downward bins. The indices L 
and H stand for the lower (E < E t h) and higher (E > Ef^) energy bins, and the threshold 
energy _E t h is chosen so that the numbers of events for the lower and higher energy bins 
are approximately equal, and in the case of the sub-GeV events, E? t h = 0.5GeV, and in the 
case of the multi-GeV events, the threshold energy is E t h = 3.2GeV, respectively, for all 
the zenith angle bins. We have put the systematic errors as follows: 


<73.sl = vpmi = 0.03, <Tfjs 2 = 0.05, crpm 2 = 0.16, a a = 0.2, (3.4) 

(jyfA = 0.005, (J 7 l2 = 0.008, (Try hi = 0.021, ct 7 h2 = 0.018, cr 7 i = 0.015, er 7 2 = 0.025. 

(3.5) 
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Figure 2. Upper panel: The allowed region in the (e ee , |e eT |) plane from the HK atmospheric 
neutrino data for a normal mass hierarchy (left panel) and for an inverted mass hierarchy (right 
panel) from the energy-rate analysis. Lower panel: The same allowed region as the upper panel 
from the two energy-bin analysis. 


All the systematic errors in (3.4) except ijg m 2 and <t 7 2 are the same as those in (2.4) in sect. 2 
and those used in Ref. [51]. crg^ = 0.16 is the uncertainty in the relative normalization 
between the multi-GeV v - v flux and it was 0.05 in (2.4). <t 7 2 = 0.025 is the uncertainty 
in the relative normalization between the upward and downward going multi-GeV /r-like 
events and it was 0.008 in the analysis of SK data [51]. The choice of these systematic errors 
(3.4) and (3.5) and the index jR = 3 has been made so that the result of our analysis on 
the mass hierarchy is close to that in Ref. [24], and we have checked that it dose not affect 
the sensitivity to NSI significantly. 

The results from the energy rate (spectrum) analysis are given by the upper (lower) 
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panel in Fig. 2. From the energy rate analysis we have |e eT /(l + e ee )| < 0.3 at 2.5erCL. On 
the other hand, from the energy spectrum analysis we get —0.1 < e ee < 0.2 and |e eT | < 0.08 
at 2.5cr (98.8%) CL for the normal hierarchy and to —0.4 < e ee < 1.2 and |e er | < 0.34 at 
2.5cr (98.8%) CL for the inverted hierarchy. 

From Fig. 2 we note two things. Firstly, the allowed regions from the energy spectrum 
analysis (the lower panel) are much smaller than those from the energy rate analysis (the 
upper panel) for both mass hierarchies. Secondly, the allowed regions (the right panel) for 
the inverted hierarchy are wider than those (the left panel) for the normal hierarchy for 
both rate and spectrum analyses. 

To understand these phenomena, we have plotted in Fig. 3 Xmuiti-GeV f° r the mu ^i- 
GeV events, which are expected to be sensitive to the matter effect and therefore to e ee , 
as a function of e ee in the case of e er = 0. In plotting the figures in Fig. 3, we have taken 
into account only the statistical errors for simplicity, and we assume that the HK detector 
could distinguish neutrinos and antineutrinos for both e-like and ^u-like events in all the 
energy ranges of the multi-GeV events, and that the detection efficiency is the same for 
both neutrinos and antineutrinos. Since the SK collaboration distinguish neutrinos and 
antineutrinos only for the multi-GeV e-like events [47], our assumption here may not be 
realistic, and the separate plots for neutrinos or for antineutrinos except for the e-like 
events should be regarded as information for theoretical consideration. The two figures 
((a) and (b)) in the top row are the results of the energy rate analysis. The two figures 
((c) and (d)) in the middle row are the results of the energy spectrum analysis with two 
energy bins for the separate neutrino or antineutrino events. The two figures ((e) and (f)) 
in the bottom row are the results of the energy spectrum analysis with two energy bins of 
neutrinos and antineutrinos combined. Comparing the figures ((a) and (b)) in the top row 
and those ((e) and (f)) in the bottom row, we see that, even if some of the data set in the 
spectrum analysis have a sensitivity to the effect of e ee , the data in the rate analysis does 
not necessarily have a sensitivity to e ee particularly for e ee > 0, for both mass hierarchies. 
While it is not clear to us why the sensitivity is lost only for e ee > 0, we have found that, 
if we try to fit the same data with the numbers of events with the wrong mass hierarchy, 
then the plot becomes left-right reversed, i.e., the sensitivity is lost only for e ee < 0. On 
the other hand, by comparing the figures ((c) and (d)) in the middle row and those ((e) 
and (f)) in the bottom row, we see that, in the case of the inverted mass hierarchy, even 
though the separate data has a sensitivity to e ee , the combined data loses a 

sensitivity. We could not explain these phenomena using the analytic expression for the 
oscillation probability, but we interpret this loss of sensitivity as a destructive phenomenon 
between neutrinos and antineutrinos in the rate analysis, and between the lower and higher 
energy bins in the spectrum analysis for the inverted mass hierarchy. 

It is expected that the HK experiment will be able to use information on the energy 
spectrum, so we believe that the allowed region in the lower panel in Fig. 2 with the energy 
spectrum analysis reflects the true HK sensitivity more than that in the upper panel does. 
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Figure 3. The behaviors of X^uiti-GeV f° r teT = 0 as a function of e ee - (a), (b): Energy rate 
analysis for NH (a) and IH (b). (c), (d): Energy spectrum analysis for NH (d) and IH (e) for the 
separate neutrino or antineutrino events, (e), (f): Energy spectrum analysis for NH (e) and IH (f) 
using only the combined numbers of events of v e + v e and In (a), (b), (c) and (d), the 

plots for the separate neutrino or antineutrino events are created based on the assumption that HK 
could distinguish neutrinos and antineutrinos. 
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3.2 The case in the presence of NSI 

Next let us discuss the case where “the experimental data” is the one obtained with 
(e ee ,e eT ) / (0,0). The analysis is the same as the one in subsect. 3.1, except that the 
“the experimental data” is produced assuming the presence of NSI, and here we perform 
only an energy spectrum analysis with two energy bins. The results are shown in Fig. 4, 
where the allowed regions at 2.5 <tCL (Ay 2 = 8.8 for 2 degrees of freedom) around the 
true points are depicted. The straight lines |e er | = 0.8 x |1 + e ee | in Fig. 4 stand for the 
approximate bound from the SK atmospheric neutrinos in Fig. 1, and we have examined 
only the points below these straight lines. As seen from Fig. 4, the errors in e ee and |e er | 
are small for |e ee | < 2 in the case of the normal hierarchy and for — 3 < e ee < 1 in the case 
of the inverted hierarchy. The errors are large otherwise, and the reason that the errors are 
large is because a sensitivity is lost due to a destructive phenomenon between neutrinos 
and antineutrinos as was discussed in subsect. 3.1. 

We note in passing that there are a couple of points in Fig. 4, where the allowed region 
has an additional isolated island. This is regarded as so-called parameter degeneracy [52- 
55] in the presence of the NSI. Since little is known about parameter degeneracy in the 
presence of the new physics and since the study of the subject is beyond the scope of this 
paper, we do not discuss parameter degeneracy here. 
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Figure 4. The allowed region at 2.5erCL around the point (e ee , \e eT \) ^ (0, 0), where 5 = arg(e eT ) = 
0 is assumed. Most of the allowed regions are connected, but those around a few points have an 
isolated island, and they are depicted with different colors. The allowed regions at e ee = 3 for the 
inverted mass hierarchy are much wider compared with other cases, so their boundary as well as 
their center are shown with dashed lines and with different colors. 
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4 Conclusions 


In this paper we have studied the constraint of the SK atmospheric neutrino data on 
the non-standard flavor-dependent interaction in neutrino propagation with the ansatz 
(1.9). From the SK atmospheric neutrino data for 4438 days, we have obtained the bound 
| ^er l/|i + fee| iS at 2 .o(jCL, while we have little constraint on c ee . 

We have also discussed the sensitivity of the future HK atmospheric neutrino experi¬ 
ment to NSI by analyses with the energy rate and with the energy spectrum. If nature is 
described by the standard oscillation scenario, then the HK atmospheric neutrino data will 
give us the bound |e eT |/|l + e ee | < 0.3 at 2.5crCL from the energy rate analysis, and from 
the energy spectrum analysis it will restrict e ee to —0.1 < e ee < 0.2 and |e er | < 0.08 at 
2.5(7 (98.8%) CL for the normal hierarchy and to —0.4 < e ee < 1.2 and |e er | < 0.34 at 2.5cr 
(98.8%) CL for the inverted hierarchy. On the other hand, if nature is described by NSI 
with the ansatz (1.9), then HK will measure the NSI parameters e ee and |e er | relatively 
well for \e ee \ < 2 in the case of the normal hierarchy and for —3 < e ee < 1 in the case of 
the inverted hierarchy. 

We have shown that it is important to use information on the energy spectrum to 
obtain strong constraint, because a sensitivity to NSI would be lost due to a destructive 
phenomena between the low and high energy events. If there is a way to distinguish 
between neutrinos and antineutrinos, as is done by the SK collaboration [47] for the e-like 
multi-GeV events, also for the multi-GeV /x-like events, then the sensitivity to NSI would 
be greatly improved, because in this case we can avoid a destructive phenomena between 
neutrinos and antineutrinos. 

While HK is expected to play an important role in measurement of 5 in the standard 
three-flavor scenario using the JPARC beam, HK has also a potential for new physics with 
atmospheric neutrinos. Search for NSI may lead to physics beyond the Standard Model, 
and the effects of NSI at HK deserves further studies. 
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